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ABSTRACT 

We report the discovery of an extremely close white dwarf plus F dwarf main-sequence 
star in a 12 hour binary identified by combining data from the RAdial Velocity Experi¬ 
ment (RAVE) survey and the Galaxy Evolution Explorer (GALEX) survey. A combination of 
spectral energy distribution fitting and optical and Hubble Space Telescope ultraviolet spec¬ 
troscopy allowed us to place fairly precise constraints on the physical parameters of the bi¬ 
nary. The system, TYC 6760-497-1, consists of a hot Te® ~ 20, OOOAT, M^d O.GMq white 
dwarf and an E8 star (Mms ^ 1.23Mq, i?MS l.SR©) seen at a low inclination (* ^ 37°). 
The system is likely the descendent of a binary that contained the E star and a ~2M0 A-type 
star that filled its Roche-lobe on the thermally pulsating asymptotic giant branch, initiating a 
common envelope phase. The E star is extremely close to Roche-lobe filling and there is likely 
to be a short phase of thermal timescale mass-transfer onto the white dwarf during which sta¬ 
ble hydrogen burning occurs. During this phase it will grow in mass by up to 20 per cent, until 
the mass ratio reaches close to unity, at which point it will appear as a standard cataclysmic 
variable star. Therefore, TYC 6760-497-1 is the first known progenitor of a super-soft source 
system, but will not undergo a supernova la explosion. Once an accurate distance to the sys¬ 
tem is determined by Gaia, we will be able to place very tight constraints on the stellar and 
binary parameters. 
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1 INTRODUCTION 

The unique properties of Type la Supemovae (SN la) as distance 
indicators, suffi ciently bright to serve as ya rdsticks on cosmolog¬ 
ical scales (e.g. [Branch & Tammanniri992h . has resulted in them 
becoming some of the most important objects in the Un iverse, hav¬ 
ing led to the discovery of i ts accelerating expansion dRiess et al.l 
ll998l : IPerlmutter et al.lll999l) . 

Although it is generally accepted that SN la are related to the 
thermonuclear ignition of a white dwarf that surpassed the Chan¬ 
drasekhar mass limit, there is no consensus yet on the pathways 
leading to the explosion. The two main formation channels are 
thought to be the “single-degenerate” channel, where a white dwarf 
accretes materi al from a non-degene rate companion yia Roche- 
lobe overflow d Whelan & IbenI ri973h during a phase of thermal 
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time-scale mass transfer known also as the super soft source (SSS) 
phase, and the “double-degenerate” channel in wh ich the explosio n 
is the result of the merger of two white dwarfs dWebbinklll984h . 
In recent years the idea of a “double-detonation” scenario has also 
been advanced. In this case a layer of helium on the surface of the 
white dwarf detonates, triggering the explosion of the underlying 
core either directly, or via a comp ressional shock wave dFink et aU 
l200ll20ld:rShen & Moorel2014l) . hence it is possiblevia this chan¬ 
nel to detonate the white dwarf while its mass is still below the 
Chandrasekhar limit. In theory these “sub-Chandrasekhar SNIa” 
are possible via both main channels. In the single degenerate chan¬ 
nel the white dwarf accretes heliu m-rich material from a dono r 
star, building up a large surface layer dTutukov & Yungelsonll99^ . 
whilst the merger of a carbon-oxygen core white dwarf with a low- 
mass helium core white dwarf can lead to heliu m ignition and hen ce 
a supernova via the double degenerate channel l lFink et aT]|2007h . 

The search for Galactic SNIa progenitors from either of the 
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channels is a difficult task. In the single degenerate case very few 
SSS sys tems sufficien tly nearby for detailed parameter studies are 
known iGreineill2000h . This is due to a combination of the short 
duration of the SSS phase and because their very soft X-ray emis¬ 
sion is easily absorbed by neutral hydrogen in the galactic plane. 
Whilst double white dwarf binaries are intrinsically faint objects 
and difficult to distinguish from sin gle white dwarfs withou t dedi¬ 
cated spectroscopic monitoring (see Napiwotzki et aljTOOl for ex¬ 
ample). However, in both cases the progenitor systems are the de- 
scendents of detached white dwarf plus F, G or early K type main- 
sequence star companions. With no ongoing accretion these sys¬ 
tems are e asy to characterise. F urthermore, they are expected to be 
numerous dHolberg et al.ll20f^ and so are optimal for population 
studies and hence testing both SN la progenitor channels. 

However, identifying white dwarfs with early type compan¬ 
ions has been extremely difficult until recently. This is due to the 
fact that the main sequence st ar completely outshines the whit e 
dwarf at optical wavelengths dRebassa-Mansergas et al.l [2012ah . 
hence a combination of optica l and ultraviolet (UV ) coverage is 
required to find these systems. iMaxted et alj d2009h attempted to 
detect systems of this type by combining data from the Galaxy 
Evolution Explorer (GALEX) survey and the Sloan Digital Sky 
Survey (SDSS), but were limited by the optical colour selection of 
their main-sequence stars. Whilst lB^urleigh et al.l d 19971) found only 
four unresolved white dwarf plus EGK main-sequence star systems 
searching for extreme-UV and soft X-ray excesses. 

We have recently begun a project combining the large 
dataset of the RAdia l Velocity Experiment (RAVE) survey 
dKordopatis et al.ir2013h in the southern hemisphere and the LAM- 
OST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope) 
survey in the northern hemisphere, with UV data from the GALEX 
survey in order to identify main-sequence F, G and early K stars 
with significant UV excesses. 

The RAVE survey spectroscopically observed more than 
400,000 bright (8< I <12 mag) southern hemisphere stars in the 
spectral region 8410-8794A (the infrared calcium triplet) with a 
resolution of R~7000. The data were used to establish the basic 
parameters of each star (effective temperatures, surface gravities 
and metallicities) as well as line-of-sight velocities. We used these 
data, in conjunction with PHOENIX synthetic spectra dHusser et al.l 
l2013h . to estimate the UV colours of these stars. We selected those 
showing UV excess flux for radial velocity follow-up observations 
to detect close binarity, as well as Hubble Space Telescope UV 
spectroscopy to confirm the presence of a white dwarf. Eull details 
of this process and the list of UV-excess objects will be published 
in a forthcoming paper (Parsons et al. in prep.) 

As well as providing a crucial test of the single degener¬ 
ate channel for SN la, these detached white dwarf plus early-type 
main-sequence star binaries also provide useful constraints for the 
common envelope phase of binary evolution, which they passed 
through when the progenitor of the white dwarf evolved off the 
main-sequence. Essentially all current observational constraints on 
the common envelope phase are based on systems containing white 
dwarf s with low-mass M dwarf companions (e.g. IZorotovic et al.l 
l20l^ . However, the handful of systems with more massive main- 
sequence star components already hint that additi onal energy may 
be ne eded to expel the envelope in these cases dZorotovic et all 
l2014h . 

In this paper we present data for the first of our UV-excess 
objects from the RAVE sample to be confirmed as a close white 
dwarf plus main-sequence binary. The object, TYC 6760-497-1, is 
a late-E star. We derive clear observational constraints on the stellar 


Table 1. Radial velocity measurements for the main-sequence star in 
TYC 6760-497-1. 


BJD 

(mid-exposure) 

Velocity 

(kms“^) 

Uncertainty 

(kms“^) 

Telescope/ 

instrument 

2456811.5672390 

30.33 

0.50 

Du Pont/Echelle 

2456827.6293689 

-48.79 

0.48 

MPG2.2/FEROS 

2456828.4841811 

53.36 

0.28 

MPG2.2/FEROS 

2456828.5729529 

-13.12 

0.28 

MPG2.2/FEROS 

2456828.5904822 

-26.22 

0.28 

MPG2.2/EEROS 

2456828.6012976 

-34.38 

0.27 

MPG2.2/EEROS 

2456828.6663049 

-61.48 

0.30 

MPG2.2/EEROS 

2456828.6775301 

-63.33 

0.30 

MPG2.2/EEROS 

2456828.7213267 

-55.85 

0.28 

MPG2.2/EEROS 

2456829.5163187 

32.99 

0.29 

MPG2.2/FEROS 

2456829.5491571 

7.39 

0.27 

MPG2.2/FEROS 

2456829.5939247 

-29.35 

0.30 

MPG2.2/EEROS 

2456829.6737526 

-62.99 

0.31 

MPG2.2/EEROS 

2456829.6855582 

-61.90 

0.29 

MPG2.2/EEROS 

2456829.6963120 

-60.80 

0.29 

MPG2.2/EEROS 

2456829.7123239 

-58.27 

0.29 

MPG2.2/EEROS 

2456829.7603942 

-33.58 

0.30 

MPG2.2/EEROS 

2456829.7718605 

-24.86 

0.29 

MPG2.2/EEROS 

2456829.7818813 

-15.99 

0.28 

MPG2.2/EEROS 

2456834.6122883 

-49.91 

0.30 

MPG2.2/EEROS 

2456834.6226973 

-52.15 

0.29 

MPG2.2/EEROS 

2456835.7338779 

-41.03 

0.29 

MPG2.2/EEROS 


and binary parameters of the systems, reconstruct its evolutionary 
history, predict its future and discuss the implications for our un¬ 
derstanding of close compact binary star evolution. 


2 OBSERVATIONS AND THEIR REDUCTION 

2.1 Optical echelle spectroscopy 

2.1.1 FEROS and the Du Pont echelle 

We obtained high resolution spectroscopy of TYC 6760-497-1 
with the echelle spectrograph (R~40,000) on the 2.5-m Du Pont 
telescope located at Las Campanas Observatory, Chile and with 
EEROS (R~48,000) on the 2.2-m Telescope at La Silla, Chile. 
EEROS covers the wavelength range from ~3500A to ~9200A. 
The observations formed part of a program to detect binarity among 
a large number of potential white dwarf plus main-sequence bina¬ 
ries. However, it was immediately clear that TYC 6760-497-1 had a 
very short period when we observed a very large velocity shift be¬ 
tween our first two measurements. Therefore, we obtained several 
spectra of this target per night throughout the rest of the observing 
run with the aim of measuring its period. 

Data obtained with the EEROS and duPont spectrographs 
were extracted and analysed with an automated pipeline developed 
to process spectra coming from different instr uments in an homoge¬ 
neous and robust manner jjordan et alj20l4 Brahm et al in prep.). 
After performing t ypical image r eductions, spectra were optimally 
extracted following lMarshI Jl989l) and calibrated in wavelength us¬ 
ing reference ThAr Lamps. Eor EEROS data, the instrumental drift 
in wavelength through the night was corrected with a secondary 
fiber observing a ThAr lamp. In the case of the duPont data, ThAr 
spectra were acquired before and after each science observation. 
Wavelength solutions were shifted to the barycenter of the solar 
system. 
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Figure 1. Obsei-ved FEROS spectrum (black line) of TYC 6760-497-1 with 
the best fit model overplotted (red line, Urot sin* = 75kms~^). We also 
plot the model spectra without any rotational broadening (grey line) to high¬ 
light the rapid rotation of the star. 


2.1.2 X-shooter 

A medium resolution spec trum (R^5,000) of TYC 6760-497-1 was 
obtained with X-shooter l lD’Odorico et alJboOhl) mounted at the 
Cassegrain focus of VLT-UT2 at Paranal on the 26th of Febru¬ 
ary 2015. X-shooter is comprised of three detectors that enable 
one to obtain simultaneous data from the UV cutoff at 0.3/im to 
the if-band at 2.4)rm. Our data consisted of three 60 second ex¬ 
posures, which were reduced using the standard pipeline release 
of the X-shooter Common Pipeline Library (CPL) recipes (version 
2.5.2). The instrumental response was removed and the spectrum 
flux calibrated by observing the spectrophotometric standard star 
CD-38 10980 and dividing it by a flux table of the same star to pro¬ 
duce the response function. The spectrum was extinction corrected 
but not corrected for telluric features. 


2.2 Ultraviolet spectroscopy 

TYC 6760-497-1 was observed with the Fluhble Space Telescope 
(HST) on the 9th of January 2015 with the Space Telescope Imag¬ 
ing Spectrograph (STIS) as part of program GO 13704. We used 
the G140L grating centered on 1425A, for one spacecraft orhit, re¬ 
sulting in a total exposure time of 2381 seconds. The data were 
processed using CALSTIS V3.4. We de-reddened the STIS spec¬ 
trum using a value of E{B — V) — 0.103, determined from our 
SED fit to the main-sequence star (see Section ldTt . 

3 RESULTS 
3.1 Optical spectra 

Radial velocities were computed from our optical echelle spec¬ 
tra using the cross-correlation technique against a binary mask 
representative of a G2-type star. The uncertainties in radial ve- 
locity were compu ted using scaling relations (for more detail see 
llordan et al.ll2014h with the signal-to-noise ratio and width of the 
cross-correlation peak, which were calibrated with Monte Carlo 
simulations. The results are listed in Table [T] We note however, 
that these relations were calibrated using slowly rotating stars and 



Frequency (cycles/day) 

Figure 2. plotted against orbital period for the fit to the radial velocity 
measurements. There is a clear minimum at a period of 0.5 days 


hence the uncertainties may he underestimated for rapidly rotating 
stars, such as the one in TYC 6760-497-1. 

We also estimated the projected rotational velocity (Vrot sin i) 
of the main-sequence star in TYC 6760-497-1 by comparing the 
obs erved FEROS spect ra against a synthetic grid of stellar spec¬ 
tra jCoelho et alj|200^ . The synthetic spectra were degraded to 
the resolution of FEROS by convolving them with a Gaussian 
with R = X/SX — 53000 and then further degraded to differ- 
ent va l ues of Wrot sin i using the limb-darkening coefficients of 
IClaretl ( l2004h . The optimal fit (which also yields a set of stellar 
parameters) was found by chi-square minimisation in 3 echelle 
order of the spectra which include the zone of the magnesium 
triplet (5000-5500A). The measured rotational broadening was 
Vrot sini = 75 ± 3kms“^. We also obtained the following stel¬ 
lar parameters: TeS = 5750 ± 200 K, log g = 3.75 ± 0.3 (in cgs 
units), [Ee/H]=—0.5 ± 0.2. Figure [T] shows the observed spectra 
with the best fitted model, as well as the same model with no ro¬ 
tational broadening. This high rotation rate means that the fitted 
stellar parameters from the spectral fit are not necessarily reliable 
since most of the weak lines that have a log g dependence are com¬ 
pletely smeared out. Moreover, the system is extremely close to 
Roche-lobe filling (see Section |4j and hence the effects of gravity 
darkening and Roche-distortion likely have an effect on the spectral 
fit. 


3.2 Orbital period 

We measured the binary period of TYC 6760-497-1 by fitting a con¬ 
stant plus sine wave to the velocity measurements over a range of 
periods and computing the of the resulting fit. The result of this 
is shown in Figure and shows a clear minimum in at a period 
of 0.5 days. We determine the ephemeris of TYC 6760-497-1 as 

BJD = 2456823.81993(40) + 0.498688(26)£', (1) 

where E is the binary phase and phase 0 occurs at the conjunction 
of the main-sequence star. 

The phase-folded radial velocity plot is shown in Figure[3 We 
find a radial velocity semi-amplitude for the main-sequence star 
of 65.0 ± 0.3kms“^and a systemic velocity of 1.4 ± 0.2kms“'^. 
The lower panel of Figure shows the residuals to the fit, which 
show variations larger than their uncertainties. This is likely due to 


© 2015 RAS, MNRAS POQ.fTIfT^ 
















4 S. G. Parsons et al. 



Orbital phose 


Figure 3. Phase-folded radial velocity plot for the main-sequence star in 
TYC 6760-497-1. The lower panel shows the residuals to the best fit. 


the main-sequence star’s large rotational broadening causing small 
systematic errors during the cross correlation process. 


4 SYSTEM PARAMETERS 

In this section we detail how we constrained the physical parame¬ 
ters of the binary. 


4.1 SED fit 

In order to get an initial estimate of the physical parameters of the 
main-sequence star in TYC 6760-497-1 and the reddening towards 
the system we fitted its spectral energy distrib ution (SED) usin g 
the virtual observatory SED analyzer iVOSA: feavo et al]|2008h . 
We convolved our flux calibrated X-shooter spectrum with a num¬ 
ber of generic narrow and broad band filters including Strdmgren 
uvby, Bessell BVRI Johnson BVR and Cousins RI filters. We 
comb ined these with archival data from the GALEX (|Martin_e^^ 
2005h . 2MASS dSkrutskie et al.l l2006h and WISE dWright et al. 
20101) catalogues. We excluded the GALEX data when fitting 
the SED since we expected the white dwarf to contribute a non- 

negligible amount of flux at these wavele ngths. _ 

We used a large grid of BTSettl d Allard et all 12012h mod¬ 
els to determine the main-sequence star’s parameters and allowed 
the extinction to vary from zero to three times the maxim um ex¬ 
pected value of Av = 0.36 dSchlaflv & Einkbeineill20ll]) . As ex¬ 
pected the fit was insensitive to the metallicity of the system and 
only mildly sensitive to the surface gravity, preferring a value of 
logg ~ 4.0. However, we found that the temperature was well 
constrained to T^s = 6400 ± 100 K and the extinction to be 
Av = 0.32, implying E{B — V) — 0.103, in good agre ement 
with what is expected from the lSchlaflv & Finkbeineil d201lh maps. 
The uncert ainties were deter mined via a Bayes analysis (for more 
details see [Bavo et al.l[200^ . These measurements are consistent 
with those from the RAVE survey, in which TYC 6760-497-1 was 
determined to have Teff = 6150 ± 100 K and logg — 3.6 ± 0.2 
dKordonatis et al.ll2013h . although these measurements are likely 
affected by the large rotational broadening, similar to our EEROS 



Wavelength (microns) 


Figure 4. The spectral energy distribution of TYC 6760-497-1 in the UV 
and optical wavelength ranges (the SED fit also included infrared data, but 
we limit the plot to this range to demonstrate the relative contributions of 
the two stars). The black lines are the observed UV HST/STIS spectram and 
optical X-shooter spectram (no telluric correction was applied). Broad band 
photometry is shown in red (their errors are too small to see on this scale), 
the best fit BTSettl main-sequence star model is shown in grey and was 
rotationally broadened to match the measured value for the main-sequence 
star in TYC 6760-497-1. We also plot in dark blue a TLUSTY/SYNSPEC 
model for a Tgu = 20400 K, log g = 8 white dwarf. The fit to the UV part 
of the spectrum is better illustrated in Figured 


fit. In Figure|4]we show part of the SED in the UV and optical wave¬ 
length ranges with the best fit model and observed spectra overplot¬ 
ted. While the GALEX NUV measurement shows little evidence of 
an excess, the FUV measurement is clearly dominated by the white 
dwarf. 


4.2 Rotational broadening constraint 


Assuming that the main-sequence star is tidally locked to the white 
dwarf, which it should be since the tid al synchroni sation timescale 
for this system is less than 0.5 Myr dZahnl [l97% . much shorter 
than the white dwarf’s cooling age (see Section ISTt . the rotational 
broadening measurement can be used to place some constraints on 
the binary parameters via 

t;rotsin2 = Kms{1 (2) 

a 

where Kms is the radial velocity semi-amplitude of the main- 
sequence star, q = Mms/AJwd, the mass ratio of the binary and 
Rms/cl is the radius of the main-sequence star scaled by the orbital 
separation. Combining this with Kepler’s third law. 


P'GMms 

47r2 V 9/ ’ 


(3) 


where P is the orbital period and G is the gravitational constant, 
we can effectively solve for q (hence Mwd) by assuming a mass 
and radius for the main-sequence star. Once this is known all the 
other binary parameters can then be determined. For example, the 
orbital inclination via 


• 3 . 
sm t 


P^Ms(i + g)^ 

2'iiGMwu 


(4) 


Ideally the fit to the EEROS spectra would yield the main- 
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Figure 5. Limits on the physical parameters of the main-sequence star. The 
measured rotational broadening, period and radial velocity amplitude ex¬ 
clude the light grey areas as at low surface gravities the star fills its Roche 
lobe, whilst high surface gravities lead to unphysical results. The horizonal 
blue lines indicate the possible temperature range from the SED fit. The 
dark grey hatched area indicates the parameter range that is consistent with 
the fit to the white dwarf’s spectrum (see Section|43). The small region that 
is consistent with all our constraints is highlighted by a red star symbol. 


sequence star’s effective temperature, surface gravity and metal- 
licity, which can t hen be used to esti mate its mass and radius via 
the Torres relation ^Torres et aljboioll and hence determine the bi¬ 
nary parameters. Unfortunately, due to its rapid rotation the FEROS 
fit does not give reliable stellar parameters. Therefore, keeping the 
metallicity fixed at the solar value, we used a grid of Teft, log g val¬ 
ues and the Torres relation to determine the range of possible binary 
parameters. The result of this is shown in Figure^ We found that at 
surface gravities below ~4.3 the binary parameters were such that 
the main-sequence star should be Roche-lobe filling. Since there 
is no evidence for this, we can exclude this region. Furthermore, 
at very high surface gravities (>4.8), a solution to the binary pa¬ 
rameters is not possible since it would require the inclination to 
be larger than 90 degrees. Hence this region can also be excluded. 
However, this still results in a large range of possible parameters, 
the white dwarf mass can be anything between 0.3Mq and O.VMq, 
the main-sequence star mass between l.OSM© and 1.26Mq and 
the inclination larger than 32 degrees. Varying the metallicity does 
not have a large effect on these limits. 

We estimated the distance to TYC 6760-497-1 using 
isochrones from the PAdova an d tRieste Stellar Evolution Code 
(PARSEC) teressan et alj|201^ . For a given temperature we used 
the isochrones to calculate the absolute magnitude of the main- 
sequence star over the full range of allowed surface gravities in 
the BVRIJHK bands and compared these to the measured val¬ 
ues to calculate the distance. This results in a range of distances 
for a given main-sequence star temperature. We take into account 
the variations in the calculated distances from the different bands 
and the effects of the unknown age of the star. The result of this is 
shown in Figure[6l While the unknown age of the star generally has 
little effect on the predicted distance, if it is very young (<30Myr, 
which is in fact ruled out by the cooling age of the white dwarf) or 
very old (>3.2Gyr) the resultant distance is somewhat larger than 
whilst on the main-sequence. 


Figure 6. Distance to the main-sequence star in TYC 6760-497-1 as a func¬ 
tion of its effective temperature. The grey region shows the permitted range. 


4.3 White dwarf spectral fit 


We fitted the HST/STIS spectrum of the white dwarf in TYC 6760- 
497-1 with a grid of white dwarf mod els computed using 
TLUSTY/SYNSPEC faubenv & Lanj|l995h . stepping in surface 
gravity from 7.0 to 9.5 in steps of 0.1 (see Figure|7]l. At each step 
the best fit to the spectrum gives the tem perature and hence th e 
mass and radius using the cooling models of iFontaine et alj ( l200lh . 
Scaling the flux then also gives an estimate of the distance. Hence 
we are able to determine the white dwarf mass, effective tempera¬ 
ture and distance as a function of the surface gravity. These rela¬ 
tions are illustrated in Figure Teff and log g are correlated pa¬ 
rameters, so increasing logf/ results in a higher best-fit Teff. At 
the same time, increasing log g implies a higher mass, and, via the 
mass-radius relation, a smaller radius. This drop in radius domi¬ 
nates over the increase in Teff in the UV flux emitted by the white 
dwarf, i.e. increasing the surface gravity leads to a lower best-fit 
distance. 

Several narrow absorption features are seen in the STIS spec¬ 
trum, primarily from Si and C as a result of the white dwarf ac¬ 
creting material from the wind of the main-seque nce star as is 
commonly seen in other close white dwarf b i naries (fTappert et al.l 
I 2 OIII; IParsons et al.ll2012l: IPvrzas et al.|[2012l: iRibeiro et alJl2013h . 


We found setting the metal abundances in the model to 0.1 times 
solar for all elements fits these well. 


4.4 Combined constraints 

We can combine the constraints from fitting the white dwarf’s spec¬ 
trum with those from the main-sequence star to determine which 
set of parameters are consistent with all our data. For a given main- 
sequence star effective temperature, we have a range of possible 
distances from the isochrone fitting (see Figure , which then 
leads to a range of possible white dwarf masses from the relations 
shown in Figure The rotational broadening measurement also 
yields a white dwarf mass for a given main-sequence star temper¬ 
ature and gravity, allowing us to determine the range over which 
these two measurements are consistent. This range is shown in Fig¬ 
ure [5] as the dark grey hatched region. We can instantly exclude a 
main-sequence star with a temperature less than ~6200 K, since a 
star this cool would be at a short distance and hence would have 
a white dwarf with a mass>O.7M0, which is only possible if the 
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Figure 7. HST/STIS spectrum of the white dwarf in TYC 6760-497-1 with 
a model fit overplotted (red line, log <; = 8.0). The dashed parts were not 
included in the ht, these include the core of the Lya line, which is contam¬ 
inated by geocoronal emission, and the red end of the spectrum, which has 
low signal-to-noise. There are also several emission lines including the quite 
strong Civ 1550A line. These are likely due to the white dwarf capturing 
material from the wind of the F star. 


system fills its Roche lobe. It also means that the system must be at 
a minimum distance of 250 pc. 

These constraints are also consistent with those from the main- 
sequence star temperature determined from the SED fitting (Tefi = 
6400 ± 100 K), which is also illustrated in Figure|^ There is there¬ 
fore a small region in which all of our measurements and fits are 
consistent, highlighted by a red star in Figure |5] The full set of 
ranges for all the parameters are given in Table |2] 


5 DISCUSSION 

Having determined the stellar and binary parameters of TYC 6760- 
497-1, we can now investigate the evolution of the system and its 
implications for models of compact binary star evolution and SN la 
formation channels. We start by reconstructing the systems evolu¬ 
tionary history. 


5.1 Constraints on common envelope evolntion 

Common envelope evolution is typically described with a 
parametrized energy equation: 

77bind (5) 


where acs represents the fraction of orbital energy that is used to 
unbind the envelo pe usually called the common envelope efficiency 
jPaczvnskill 19761) . 

The binding energy of the envelope is often assumed to be 
equal to the gravitational energy of the envelope: 


77bind 


— i7gr — 


\Ri 


( 6 ) 


where Mi, Mi,e and Ri are the total mass, envelope mass and 
radius of the primary star, and A is a binding energy parameter. Al¬ 
though very often ignored, the binding energy parameter strongly 
depends on the mass and evolutionary state of the white dwarf pro¬ 
genitor when filling its Roche-lobe. This is especially true if the 



Figure 8. Constraints on the white dwarf’s parameters as a function of its 
surface gravity. 


Table 2. Physical and binary parameters of TYC6760-497-1. 


Parameter 

Value 

RA 

15:02:22.4896 

Dec 

-29:41:15.666 

B 

11.578 

V 

11.245 

R 

11.030 

J 

10.167 

H 

9.918 

K 

9.853 

Distance (pc) 

250-320 

Orbital period (days) 

0.498688(26) 

Ams (kms“l) 

65.0 ±0.3 

Vrot sin i (km s“ ^) 

75.0 ±3.0 

Inclination (°) 

33^3 

Separation (Rq) 

3.20-3.28 

White dwarf mass (Mq) 

0.52-0.67 

White dwarf temperature (K) 

19,500-21,000 

Main-sequence star spectral type 

F8 

Main-sequence star temperature (K) 

6300-6500 

Main-sequence star surface gravity 

4.31^.48 

Main-sequence star mass (Mq) 

1.22-1.25 

Main-sequence star radius (R©) 

1.18-1.40 


recombination energy Urec available within the envelope supports 
the ejection process. Therefore, a more general form for the binding 
energy equation is: 



Gm , 

— - rdm + Qrec 

rim) 



(7) 


where ctrec is the fraction of recombination energy that contributes 
to the ejection process. It is of outstanding importance for our 
understanding of compact binary star evolution to observationally 
constrain the values of both common envelope efficiencies and to 
investigate possible dependencies on the binary parameters. 

With its period of less than half a day TYC 6760-497-1 is the 
first short orbital period (Porb < 1 day) post common envelope bi¬ 
nary (PCEB) with a massive secondary star (spectral type earlier 
than K). The two systems that have been previously discovered are 
IK Peg and KOI-3278 with much longer orbital per iods of 21.722 
and 88.18 days respectively llWonnacott et al Jfl 99.31 : iKruse & Agoll 
l2014h . These two systems would have formed through common 
envelope evolution only if the common envelope efficiency acs 
has been larger than it seems to be required to understand PCEBs 
with M d warf secondaries, where qqe ^ 0.25 — 0.5 seems to 
work best dZorotovic et al.ll201(ll : [Rebassa-Mansergas et alj|2012d : 
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Figure 9. The mass of the white dwarf after the common envelope phase 
as a function of its initial progenitor mass. The resultant white dwarf mass 
depends heavily upon the evolutionary stage at which the common envelope 
started. We show the results for a progenitor on the first giant branch (light 
grey), on the early AGB (dashed line) and on the TP-AGB (dark grey). The 
horizontal lines correspond to the range of white dwarf masses in agreement 
with our observations. 

iToonen & Neleman3l2013l : ICamacho et alj|2014h and/or if in addi¬ 
tion to a high fraction of the released orbital energy also recom¬ 
bination energy contrib uted to some degree to the e jection process 
expelling the envelope dZorotovic et al.ll20l^l2014h . 

To investigate whether TYC 6760-497-1 confirms or disproves 
the trend of larger efficiencies being re quired to unders t and PC EBs 
with massive secondaries we follow IZorotovic et alJ i2010l) and 
reconstruct the evolutionary h istory of the syste m using the bi¬ 
nary star evolution code BSE jHurlev et alj|2003) and the param¬ 
eter constraints derived in the previous sections, i.e. Mwd = 
0.52 - 0.67Mq, Mms = 1-22 - 1.25Mq, Poa = 0.498688 days 
and Teff.wD = 19500 — 21000K. The current cooling age of the 
white dwarf is only ~ 0.1 Gyr dpontaine et al.l200ll) hence the cur¬ 
rent orbital period is virtually identical to the period the system had 
at the end of the common envelope (even if efficient magnetic brak- 
in g is assumed). Reconstr ucting the common envelope evolution as 
in izorotovic et al.l ( l2014h we find a large range of possible values 
for the common envelope efficiencies for TYC 6760-497-1 and dif¬ 
ferent possible evolutionary stages at which the progenitor of the 
white dwarf could have filled its Roche-lobe. 

Figure shows the solutions with acE < 1 in a final ver¬ 
sus initial mass plot for the primary if contributions from re¬ 
combination energy are ignored (i.e. Qrec = 0.0). Solutions ex¬ 
ist for massive progenitors on the first giant branch (FGB), on 
the early asymptotic giant branch (AGB) and on the thermally 
pulsating-AGB (TP-AGB). The more massive the progenitor, the 
larger the value of the common envelope efficiency must have 
been and the younger is the system today. The possible ranges 
are Q = 0.90 — 1.0 for a 3.40 — 4.13Mq progenitor that fills 
its Roche-lobe on the FGB (becoming first a naked helium star 
then eventually a carbon-oxygen white dwarf, Ikturlev et alj|200^ 
at an age of 0.17 — 0.27 Gyr; a = 0.08 — 0.89 if the progenitor 
reached the early AGB, after 0.48 — 2.49 Gyr and with an ini¬ 
tial mass of 1.61 — 2.99M©; and finally acE = 0.05 — 0.20 for 
a 1.37 — 2.99M0progenitor that filled its Roche-lobe on the TP- 
AGB at an age of 0.48 — 4.15 Gyr, making it the oldest option. In 
the latter case we get the largest ranges as the core mass on the 


TP-AGB reaches the current white dwarf mass for a large range 
of initial masses. We therefore consider this the most likely sce¬ 
nario. This is further supported by our measurement of the radius 
of the secondary. Assuming solar metallicity a 1.235M© main se¬ 
quence star needs ~ 1.8 Gyr to expand from its ZAMS radius to 
the current radius of the secondary of 1.35R0. We therefore con¬ 
clude that the most consistent scenario for the evolutionary history 
of TYC 6760-497-1 is that the progenitor of the white dwarf was 
of a relatively low mass (~ 2 M 0 ) and filled its Roche-lobe on 
the TP-AGB, which implies that the common envelope efficiency 
must have been small acE = 0.05 — 0.20. This is much smaller 
than values obtained for the long orbital period systems with sec¬ 
ondary stars of similar mass (IK Peg and KOT3278). In addition, 
contributions from recombination energy are not required to un¬ 
derstand the existence of TYC 6760-497-1. If recombination en¬ 
ergy contributed to expelling the envelope of the progenitor of the 
white dwarf, the fraction of orbital energy lost during the common 
envelope is reduced even further while the possible initial masses 
would remain virtually identical. It therefore seems that a simple 
dependence of the total com mon envelope efficiency on the sec¬ 
ondary mass as speculated by IZorotovic et alj i2014h remains an 
incomplete prescription for common envelope evolution and other 
parameters such as perhaps the evolutionary state of the progenitor 
when it fills its Roche-lobe may play an important role. Character¬ 
izing more PCEBs with massive secondaries is therefore urgently 
required to progress with our understanding of compact binary evo¬ 
lution which is directly related to our understanding of SN la pro¬ 
genitor systems. 

5.2 The future of TYC 6760-497-1 

While it is clear that PCEBs with M dwarf secondaries evolve into 
cataclysmic variables (CVs) as long as the mass transfer will be 
dynamically stable, the future of PCEBs with more massive sec¬ 
ondary stars is more uncertain and potentially very interesting as 
the total mass of the binary usually exceeds the Chandrasehkar 
limit. Indeed, PCEBs with massive (G or F type) secondaries are 
the progenitors of SN la explosions for both the double and the sin¬ 
gle degenerate channel. If the separation of a given PCEB is large it 
potentially survives a second common envelope forming a double 
degenerate system which may then merge and produce a SN la. If 
on the other hand, the system is close it might start thermal time 
scale mass transfer and reach the rates required for stable nuclear 
burning on the surface of the white dwarf allowing its mass to in¬ 
crease and potentially reach the Chandrasehkar limit. 

With its short orbital period and given that the system is rather 
young, TYC 6760-497-1 will clearly start mass transfer before the 
secondary evolves off the main-sequence. It is also clear that the 
system will not undergo dynamical time scale mass transfer as 
the critic al mass ratio f or dynamically unstable mass transfer is 
gcr ^ 2.5 iOe et al.l2013h for donor stars with masses and radii sim¬ 
ilar to the Sun. It is less clear, however, if the mass transfer will be 
thermally stable (i.e. stable against thermal time scale mass trans¬ 
fer). The critical mass ratio for marginal stability against thermal 
time scale mass transfer in a semi-detached binary star can be ob¬ 
tained by equating the mass-radius exponent of the Roche-lobe and 
the star, i.e. 

_ dln(i?Ms) _ dln(i?L) 

^ dln(MMs) ^ dln(MMs)’ 

(see also e.g. Ide Koo]||l992l) . where the Roche-lobe radius (Rl) is 
a function of the mass ratio and the binary separation. Despite the 
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Figure 10. Critical mass ratio for thermal time scale mass transfer as a 
function of the main-sequence star’s mass for conservative (grey) and non¬ 
conservative (black) mass transfer. The position of TYC 6760-497-1 is in¬ 
dicated in red. For completeness we also provide the limits for dynamical 
time scale mass transfer (dashed lines) where we used the fits to d etailed 
calculations of the adiabatic mass-radius exponent jHiellmindl 19891) . 


secondary in TYC 6760-497-1 being slightly evolved we can get 
a first hint about the future of the system by calculatin g the limit 
implie d by Eq. 8 using the ZAMS M-R relation from iTout et aP 
( Il996l) . For conservative mass transfer we get a critical mass ra¬ 
tio for thermal time scale mass transfer of t/cr = 1.4 — 1.48 for 
TYC6760-497-1. Thus in the case of conservative mass transfer 
the system will certainly experience thermal time scale mass trans¬ 
fer. 

However, the assumption of conservative mass transfer is not 
necessarily correct. As long as the mass transfer rate stays helow 
the critical value for stable hydrogen burning on the white dwarf, 
mass transfer will most likely not be conservative. Instead, nova 
eruptions will occur, the white dwarf mass will remain nearly con¬ 
stant and angular momentum will be taken away from the system 
hy the expelled material. In Figure [TO] we show the critical values 
for the mass ratio as a function of secondary mass for hoth con¬ 
servative and non-conservative mass transfer. In the latter case we 
assume the mass expelled during nova eruptions to carry the spe¬ 
cific angular momentum of the white dwarf. The white dwarf mass 
is assumed to he constant and the mass (and angular momentum) 
loss to he continuous which has heen shown to be equivalent to a 
discontinuous sequence of nova cyc les in terms of the secular evo¬ 
lution of CVs dSchenke r et al.l 1 9981). T hese assumptions are typical 
assumptions for CVs (e.g. iRitterl 19881) . The position of TYC 6760- 
497-1 is very close to the limit for thermal time scale mass transfer 
in the case of non-conservative mass transfer so the system may 
indeed he a progenitor of a super soft source. 

In order to investigate the future of TYC 6760-497-1 in 
more detail we pe rformed dedicated simulations using MESA 
JPaxton et alJl 201 11) . We assume the initial masses of the two stars 
to he Mms = 1.235Mq and Mwd = 0.63Mq and secondary 
radius of l.dSR©, which corresponds to the mean values of the 
ranges in agreement with our observations. The starting model was 
obtained by evolving the secondary until it reaches the required ra¬ 
dius (after 1.8 Gyr) and the white dwarf is approximated as a point 
mass. We assume mass transfer to be nearly conservative (90 per 
cent of the transferred mass remains on the white dwarf) if the mass 



0.0 2.0 4.0 6.0 8.0 10.0 12.0 

Porb (hr) 


Figure 11. The predicted future evolution of the mass transfer rate (top 
panel), the white dwarf mass (dotted line, bottom panel), and the secondary 
star mass (solid line, bottom panel) as a function of orbital period. The hor¬ 
izontal red line in the top panel represents the mass transfer rate required 
for stable hydrogen burning for the respective mass of the white dwarf (i.e. 
when the mass transfer rate is above this line stable hydrogen burning oc¬ 
curs). At Porb ~ 6.5h TYC6760-497-1 reaches this critical mass transfer 
rate. As a consequence mass transfer becomes conservative and the white 
dwarf grows in mass. At Porb ~ 5 h the white dwarf mass has reached 
O.79M0and the mass transfer rate drops below the rate required for stable 
burning. At this moment TYC6760-497-1 becomes a “normal” CV with 
non-conservative stable mass transfer driven by angular momentum loss 
only and a constant white dwarf mass. 


transfer is large enough to generate stable hydrogen burning on the 
surface of the white dwarf as we run into numerical problems in 
the fully conservative case. For mass transfer rates below the limit 
for stable burning we assume that all the accreted mass is expelled 
during nova eruptions leaving the system with the specific angular 
momentum of the white dwarf. Finally, for angular momentum loss 
due to magnetic brakin g we assume the standard prescription from 
iRappanort et al.l j 1983b . 

Based on these assumptions, the future of TYC 6760-497-1 
can be estimated. The black line in the top panel of FigurefTTIshows 
the predicted mass transfer rate of TYC 6760-497-1 as a function 
of orbital period. At an orbital period of 11.9 h (in less than 5 Myr) 
mass transfer will start. The mass transfer rate will quickly increase 
as the system is thermally unstable, but does not reach the limit for 
stable hydrogen burning until the system has an orbital period of 
6.2h (13.1 Myr from now). At this point, mass transfer becomes 
conservative (we assume that 90 per cent of the burned material 
remains on the white dwarf). As a consequence, the mass transfer 
rate increases significantly and the white dwarf mass grows (bot¬ 
tom panel). When the white dwarf mass reaches 0.79M© (14.1 Myr 
from now) the mass ratio will be close to unity (secondary mass of 
O. 8 M 0 at an orbital period of 5.1 h) and the system becomes stable 
against (conservative) thermal time scale mass transfer. Therefore 
the accretion rate drops below the value required for stable hydro¬ 
gen burning. From this moment on, TYC 6760-497-1 will behave 
as expected for normal CVs: after 105 Myr it will enter the period 
gap and restart mass transfer at a much lower rate at an orbital pe¬ 
riod of 2.1 h (1.5 Gyr from now). Finally, the system will reach the 
orbital period minimum (~1.1 h) 4.8 Gyr from now. 

While the predicted final evolution of the period and mass 
transfer rate of TYC 6760-497-1 is identical to those of CVs de- 
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Figure 12. Surface abundance ratios (C/N, dashed line and C 12 /C 13 , solid 
line) for the main-sequence star in TYC 6760-497-1 as a function of orbital 
period. Once a deep convective envelope forms (Porb ~5h) both abun¬ 
dance ratios substantially decrease as carbon depleted material is dredged 
up from the inner regions of the star where CNO burning has taken place. In 
contrast, the abundance ratios of a CV born with a low-mass donor (which 
never began CNO burning) will remain largely fixed, even when the star be¬ 
comes fully convective. We show the abundance ratios for such a “normal” 
CV (C/N, black square, C 12 /C 13 , black diamond). 


scending from PCEBs with less massive secondaries, its past as a 
PCEB containing a more massive secondary star will be imprinted 
in the relative abundances of carbon and nitrogen of the accreted 
material as soon as the outer convection zone reaches regions con¬ 
taining CNO processed material. Decreased values for C/N and 
C 12 /C 13 in CVs descending fro m thermal time scale rn ass transfer 
systems have been predicted by ISchenker et alj ( I2OO2I) and found 
in CVs based on the presence of C/N line ratios measured with 
HST by lOansicke et alJ ( l2003l) . In Eigure {T2\ we predict the C/N 
and C 12 /C 13 abundance ratios for TYC 6760-497-1 as a function 
of orbital period. As soon as the secondary star starts to develop 
a deep convective envelope both surface abundance ratios decrease 
significantly. Most dramatically, C/N decreases by several orders of 
magnitude. TYC 6760-497-1 is therefore li kely a progenitor of the 
CVs with decreased C/N line ratios found bv iGansicke et alj 1 I2OO3I) 
which are frequently called/a//ed SNIa. 

A final note concerning the future of TYC 6760-497-1 con¬ 
cerns the expected white dwarf mass when the systems becomes a 
CV. The me an white dwarf mass i n CVs derived from observations 
(~ 0 . 8 Mq . rZorotovic et alJl2oTlh is significantly higher than pre¬ 
dicted by binary population models of CVs and significantly larger 
than the mean value observed in CV prog enitors containing low 
mass secondaries (see IZorotovic et al]|201ll . for details). The value 
predicted by our simulations for TYC 6760-497-1, however, is very 
close to the mean white dwarf mass of the observed CV sample. 
One could therefore speculate that if the number of CVs descend¬ 
ing from PCEBs with massive secondaries has been underestimated 
previously, the problem with the white dwarf masses in CVs could 
be solved. Assuming that a large number of CVs are descendants 
from PCEBS with massive secondary stars such as TYC 6760-497- 
1 indeed predicts an increased number of CVs with white dwarf 
masses O. 8 M 0 . However, the large number of CVs with evolved 
secondaries predicted in this scenari o violates the gene ral expla¬ 
nation for the orbital period gap (see lwiinen et al.ll2015l . for more 
details). We therefore believe that the fraction of CVs descending 


from systems with initially massive secondaries does not signifi¬ 
cantly exceed the ~ 1 0 per cent of CVs wi th decreased C/N line 
ratios as measured by iGansicke et alj l l2003h and the white dwarf 
mass problem remains to be solved. 

We emphasize that the above evolutionary scenario for 
TYC 6760-497-1 is based on several assumptions such as the 
strength of magnetic braking, the angular momentum taken away 
from the system during nova eruptions, and the fraction of mass 
that is lost during thermal time scale mass transfer. However, our 
simulation shows that it is likely that TYC 6760-497-1 is the first 
known progenitor of a super-soft source. 


6 CONCLUSIONS 

We have identified TYC 6760-497-1 as a close binary consisting 
of a white dwarf and an E8 star that is extremely close to filling 
its Roche lobe. Assuming that the F star is tidally locked, we con¬ 
strain the masses of the two stars to Mwd = 0.52 — O.67M0, 
Mms = 1-22 — I. 25 M 0 and the radius of the F star to Rms ~ 
1.18 - 1.4OR0. The white dwarf is still hot (19, 500 - 21, OOOK) 
and so the system only emerged from the common envelope 0.1 Gyr 
ago. Its progenitor was likely on the TP-AGB when the common 
envelope began. However, this means that the common envelope 
efficiency was quite low, in contrast to other systems containing 
white dwarfs with early type companions, implying that the effi¬ 
ciency depends upon more than just the mass of the companion 
star. 

TYC 6760-497-1 will become a semi-detached system in less 
than 5 Myr and will undergo a short phase 1 Myr) of thermal 
timescale mass-transfer in 13.1 Myr resulting in a roughly 20 per 
cent increase in the white dwarf’s mass, still well short of the Chan¬ 
drasekhar limit, after which it will become a standard cataclysmic 
variable system. 

Although we have placed some constraints on the binary pa¬ 
rameters, these will be greatly improved with an accurate distance 
measurement from Gaia which should reach precision of ~ 0.3% 
at the expected distance of TYC 6760-497-10. This will dramati¬ 
cally reduce the reliance on evolutionary models in our measure¬ 
ments allowing us to more rigorously test the evolutionary scenar¬ 
ios of this (and similar) systems. 
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